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1. Introduction 

Since its first appearance [I], the Jaynes-Cummings (JC) model has been extensively 
used to study radiation-matter interaction [2] in contexts like cavity quantum 
electrodynamics (CQED) [3j HJ [5] and ion traps [6]. In a CQED scenario [3], for 
instance, the model describes the one-photon interaction between a two-level atom and 
a quantized normal mode of the electromagnetic cavity field, while in ion traps [6] it 
describes the interaction between a two-level atom and a vibrational mode of its center 
of mass. In both contexts, the quantum effects predicted by the model, such as Rabi 
oscillations and collapses and revivals of the atomic inversion operators, have been 
experimentally observed. 

Concentrating on CQED, it is of central importance to describe cavity losses, due 
for example to dissipative effects in the cavity mirrors. In Ref. [Tj we addressed this 
problem and provided a microscopic derivation for the master equation for the JC model 
neglecting spontaneous emission and supposing that the cavity was coupled to a bosonic 
reservoir. The derived master equation was then compared to the phenomenological 
master equation used in the literature [3] and some physical examples were given to 
show that the two models describe very different decay mechanisms. Indeed, while in 
the phenomenological model it is the cavity only which directly decays and the atom 
loses energy through its coupling with the cavity mode, in our model it is the atom- 
cavity system as a whole which leaks, although only the cavity is coupled to an external 
environment. 

Our approach in Ref. [7] relied on the rotating wave approximation (RWA) 
performed on the dissipator, as done in the usual formalism of master equations [8]. 
In fact it is possible to perform RWA in our case only when the vacuum Rabi frequency 
is much larger than the largest decay rate involved, i.e., in the strong coupling regime. 
Since the phenomenological model is claimed to be valid also in regimes wherein the 
vacuum Rabi frequency and the cavity decay rate are of the same order, the comparison 
between it and the microscopic model will be complete only when a microscopic 
derivation is given which takes into account in the dissipator all the terms oscillating at 
frequencies of the order of the vacuum Rabi frequency. 

This is what we are going to do in this paper. Our present investigation, based on 
a deeper analysis of the terms neglected in the so-called dressed-state approximation 
performed on the phenomenological model [3j, shows that the microscopic and the 
phenomenological model can be much closer than one might think in the context of 
our previous RWA model. In particular we will see that the two models may coincide if 
the spectrum of the environment is flat, even for small nonzero temperatures. 

The paper is structured as follows. In Sec. 2 we review the JC model, the usual 
phenomenological way to describe cavity losses and our microscopic model derived in 
Ref. [7], discussing its validity and limitations. In Sec. 3 we derive the microscopic model 
beyond RWA and in Sec. 4 the latter is compared to the phenomenological model, 
singling out the conditions under which the two models coincide, both exactly and 
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approximately. Finally in Sec. 5 some conclusive remarks are given. 



2. The Jaynes-Cummings model: the phenomenological description of 
cavity losses 

In Ref. [7] we addressed the problem of the description of the dissipation and decoherence 
processes in a two-level atom-cavity system, neglecting spontaneous emission, and taking 
into account the coupling between the cavity and an external bosonic environment. 
Denoting by \g) and |e) the atomic ground and excited states respectively, and calling 
u the atomic Bohr frequency, the atom-cavity interaction, in the RWA and in units of 
h, at resonance, is described by the Jaynes-Cummings Hamiltonian [1]: 

Hjc = ~Y a z + oo a) a + Vl (aa + + a) o _) . (1) 

Here a) (a) denotes the creation (annihilation) operator of the electromagnetic mode, 
and the atomic degrees of freedom are described by er_ = \g) (e|, a + = \e) (g\, and 
o~ z = | e) (e | — |g)(g|.The eigenstates and eigenvalues of Hjc are [2]: 

\E N>± ) = ^=(\N,g)±\N-l,e)), 

E N>± = (n - J ) lo q ± Qy/N, (2) 



2, 

except for the ground state which is given by: 

\E Q ) = \0,g), £o = ~y, (3) 

where \N,i) = \N)\i), with i = e,g, indicates the tensor product of the Fock state \N), 
with a^alN) = N\N), and the electronic state 

In the usual treatment of the problem of cavity losses one derives microscopically the 
master equation for the cavity only, when it is coupled to an environment at temperature 
T [9] and then assumes that the atom inside the cavity causes only a change in the 
Hamiltonian governing the unitary part of the dynamics [3j. In this way the following 
phenomenological master equation for the density operator of the atom-cavity system 
is assumed to be valid: 



p = - i [Hjc, p] + 7 [n(wb) + 1] 



apo) — - (a)ap + pa) a) 



a) pa — - (aa) p + paa^) 



(4) 



where n(uo) is the average number of quanta of the reservoir in the quantized mode 
of frequency Uq, and 7 is the rate of loss of cavity photons. For a reservoir at zero 
temperature, Eq. PJ becomes: 

p = -i [Hjc, p] + 7 ( apa ] - ^a f ap - ^pa^aj . (5) 
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Equations (jl]) and §5$) are the phenomenological models that are very often used in 
the literature on cavity QED pH [EH M US M UM H ■ 

In Ref. [7] we investigated the legitimacy of Eq. (TjO). By applying the general 
formalism from the theory of open quantum systems [8] , we microscopically derived the 
master equation for the system under scrutiny assuming from the very beginning that 
the Hamiltonian of the system was the JC model. 

Stated another way, the general theory of master equations in the Born-Markov 
approximation claims that the decay channels are described by jump operators which 
are transition operators between eigenstates of the Hamiltonian of the system, with 
coefficients depending on the details of the coupling between the system and the 
environment. This feature of the theory is quite general and holds even when the 
system is bipartite and only one of its two parts is coupled to an external environment. 

By applying this formalism to the JC model, we found [7] that, although cavity 
only is coupled to an external environment, it is instead the atom-cavity system as a 
whole which decays, the irreversible transitions giving rise to dissipation and decoherence 
through jumps between dressed states. We then investigated the time evolution of the 
atom-cavity system predicted by the phenomenological model and by our microscopic 
model to clarify the different physical mechanisms underlying the two different models. 

We stress that this approach includes the so-called rotating wave approximation 
(RWA) 3 on the dissipator of the master equation and is at the origin of its limited 
range of applicability expressible as [7J: 

2^ > 7max- (6) 

In other words the RWA in Ref. [7] is done over timescales of the order of the inverse 
of the Rabi frequency 2Q of the atom-cavity system, and the microscopic master equation 
is valid only when the decay occurs on a timescale which is much longer than the Rabi 
period. On the other hand the phenomenological model described by Eqs. ()!]) and (jSJ) 
uses a dissipator on which the RWA has been done over frequencies which are of the 
order of Uq, which is much larger than 2Q: therefore the phenomenological model is 
claimed to be valid in a much larger range of values of the decay rate, since in this 
model one must have 

> 7- (7) 

In the following we investigate how the master equation for the JC model is changed 
when its microscopic derivation does not neglect those terms in the dissipator oscillating 
at frequencies of the order of the Rabi frequency. 

| One should avoid confusion between the RWA performed on the Hamiltonian in the JC model and 
the RWA performed on the dissipator. In the following we will always refer to RWA as to the operation 
of neglecting the oscillating terms in the dissipator in the Interaction Picture when deriving the master 
equation. 
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3. The microscopic master equation for the JC model beyond the RWA 

In Ref. [7] we have assumed the validity of the RWA on timescales of the order of the 
period of the Rabi oscillations. Now we instead take into account the terms in the 
dissipator oscillating at frequencies proportional to Q, in order to obtain a microscopic 
model suitable for a comparison with the phenomenological model in all damping 
regimes. 

What we are going to do now is to skip the step corresponding to RWA in the 
general formalism to get a master equation under the Born and Markov approximations 
and to see what we can say about the master equation for the JC model starting from 
the following model for the total closed system: 

Hs = Hjc, H e = ^ uj k blb k , 

k 

H int = (a + a f ) 9k (b k + bfy , (8) 

k 

which describes an atom-cavity system wherein the cavity is coupled to a bosonic 
environment through the interaction hamiltonian H mt . This model is equivalent to 
that used in Ref. [9], since the additional terms rapidly oscillating at frequencies of the 
order of ujq are washed out in the total and partial RWAs performed in the following. 

From the interaction Hamiltonian H int one gets the jump operators, according to 
the recipe given in Ref. [8]: 

A (E N ^i—E N ^ m ) = \E Nt „i) (EN, m \(a + o^)\En',i) (En',i\ 

.(yW+l + lmVN) \E N>m ) (E N+1>1 \, (9) 

for N > 1 and 

A(E h± -E ) = ^=\E )(E h± \, (10) 

for N = 0. In Eq. (jHJ) we indicate the states |-E>jv,±) by |-£7jv,±i) an d the energy eigenvalues 
En,± by £?jv,±i- Accordingly / and m take the values ±1. So all the possible jump 
operators are transition operators describing jumps between eigenstates (dressed states) 
of Hjc differing for one excitation only. 

We will use the following notation for the jump operators: 

Amm = A(En+i,i — EN,m) (11) 
for N >1 and 

A 0l = A{E hl - E ) (12) 

for N = 0. 

Before explicitly inserting these operators in the master equation, let us examine 
the structure of the equation for the reduced density operator of the atom-cavity system, 
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after the Born and Markov approximations. The equation can be cast in this form: 

j t p(t)= E {e^'-^rH[AH P (t)At(^-At(^)AMp(t)]+h.c. 



u,u'>0 



+ E {e* (a,_w,)t r(-w) [A\u)p{t)A{u')-A{u')A\u)p{t)] +h.c. 

w,a/>0 

+ E {e- i( "' +w)t rV) [A(u;)p(t)A(u/) - A(w')vl(w)p(t)] + h.c.} 
w>o,a/>o 

+ E {e^' + ^r(-^)[At( w )p(t)At(cu , )-At(^)^ t (^)pW]+h.cj,(13) 

w>0,u/>0 

where we have explicitly written the minus sign for each negative Bohr frequency. The 
coefficients T(u) are complex functions of the particular Bohr frequency involved and 
are equal to: 

r+oo 

I»= / dre^(E(r)E(0)), (14) 
Jo 

with E{r) = J2 k 9k{b k e- iuJkT + b{e iuJkT ). 

All the positive Bohr frequencies involved in the master equation can be written in 
the following way: 

u Nlm = E N+l: i - E N>m = u + (ly/N + 1 - mVN^j Q, (15) 

from which, remembering that ujq ^> Q, it follows that the first two rows in the master 
equation are composed of terms which are either stationary or oscillating at frequencies 
proportional to the Rabi frequency, since: 

uJNim ~ un'm = (ly/N + 1 - m\fN - k^N' + 1 + nVW'^j Q, (16) 

while the last two rows are composed of terms oscillating at much larger frequencies, 
since: 



,kn 



2to + (l\/N + 1 - my/N + ky/N' + 1 - ny/W^j Vt. (17) 

Now we will neglect only the last two rows in Eq. (TlBl) . i.e., all the rapidly oscillating 
terms, whereas in Ref. [7] we neglected all the time-dependent terms, including the slow 
ones: we will call Quasi-RWA Master Equation the resulting equation. 

Neglecting these terms in Eq. ffTSl and going back to the Schrodinger picture, we 
then get: 

j t p(t)= -i[H JC ,p(t)} 

+ J2 {m[A(u)p(t)A\co')-A^u')A(u)p(t)]+h.c^ 

uj,uj'>0 

+ E {n-^)[AH^P(t)A^)-A(uj')A^uj) P (t)]+h.c.}. (18) 
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Inserting the explicit form of the jump operators for our model, we finally obtain: 
j t p(t)= -i[H JC ,p(t)} 



km=±l 

+ ^ ^ 0m (AomP(t)AH, kl - A^ N , kl AQ m p(t) 

N'>l;k,l,m=±l 

+ ^2 ^Nmn (A-Nm,nP{t)A^ k — A\ k AN mn p(t) 

N>l;k,m,n=±l 

+ ^ ] ^Nmn {^ANrnnP{t)A^ N i kl — Ajy; kl A]\f mn p(t) 

N,N'>l;k,l,m,n=±l 

+ f m(4mP(*)Afe-^0fc4mP(*)) 
fc,m=±l 

+ Yl ^Om {^OmPi^N'kl - A N , kl Al m p(t)) 

N'>l;k,l,rn=±l 

+ Yl ^ Nmn ( A iVmnP(*)^0fe - A 0kA ] Nmn p{t)^j 

^ f jVmn (A ] Nmn p(t)A N , k i - A Nlk iA ] Nmn p(t) 



N>l;k,m,n=±l 

+ 



N,N'>l;k,l,m,n=±l 

+ h.c, (19) 

where r 0m = T(Ei jTn — E ), T 0m = T(E — E^ m ), r mm = T(E N+ ij — E N)Tn ) and 
^Mm = r(£ , 7v,m — are complex coefficients depending on the thermal spectral 

density of the environment with which the atom-cavity system interacts, according to: 

T(oj) = ttJ(uj) + iV [ + °° dw'-^-„ (20) 

J-oo U-U' 

where the thermal spectral density J{oo) = f_°° e~ lu)T (E(t)E(O)) is the Fourier 
transform of the environment correlation function. This quantity is in general equal 

{(n(u) + 1) J Q (u) , to > 0, 
(21) 
n(\u\)J (\u\) , u<0, 

where Jq{uj) is the zero-temperature spectral density, defined for nonnegative 
frequencies, depending only on the density of modes in the reservoir and on the 

§ Eqs. ([20]) and (j!?Tj) can be obtained by generalizing the calculations given for the two-level atom 
master equation in Refs. [8l 116). 



Cavity losses for the dissipative Jaynes-Cummings Hamiltonian beyond RWA 



8 



distribution of the system-reservoir coupling constants (i.e., the g k coefficients in the 
model given in Eq. (jSJ)), and n(u) = 1/ [exp (hw/kBT) — 1] is the average number of 
photons in a mode of frequency u, for a general temperature T. The real part of T(u) is 
equal to half the decay rate relative to the Bohr frequency u, while the imaginary part 
is equal to the corresponding Lamb shift [8]. 

Equation ffl9|) is the most general equation one can give in accordance to the 
Born-Markov approximation and a RWA performed on a timescale of the order of uJq 1 - 
Anyway, Eq. ( fl9l) is not in Lindblad form and in general it might violate the complete 
positivity requirement for a master equation [8]: the possibility of casting the Quasi- 
RWA master equation in a form suitable for physical applications, i.e. Lindblad form, 
strongly depends on the specific form of the spectral density of the environment and 
must be evaluated in specific cases. We will see an important example in the rest of the 
paper. 



4. Comparison: conditions for the validity of the phenomenological model 

In order to compare the phenomenological model to the Quasi-RWA microscopic model 
given by Eq. (fl9|) . it is convenient to write Eq. in the dressed states basis. In 
the atom-cavity strong coupling regime, this procedure, combined with a secular 
approximation, leads to the so called dressed state approximation P, QUI fTT] . which 
in Ref. [7] we addressed as the main reason why the phenomenological model is so 
successful in describing the experimental situations. 

To this aim let us write the cavity mode annihilation operator with respect to the 
dressed states basis: 

oo oo 

a = ^Viv+iiiv) (jv+i|<g>i at = J2 A ° m + Yl 2 Am ™ (22) 

N=0 m=±l N=l l,m=±l 

obtained by means of Eqs. (J2J) and ([3]) and of the relation I at = (\g)(g\ + | e) (e | ) . 

Substituting Eq. (122]) and its hermitian conjugate into Eq. (j3J), one finally gets the 
phenomenological master equation in the dressed states basis: 

j t P(t)= -i[Hj C ,p(t)} 

+ 7(»N + 1) ( A °rnP(t)Al k -^{Al k A 0m ,p(t) 
_km=±l ^ 

+ ( A 0mP(t)A^ kl - ~ {A^Ao™,^)}) 

JV'>l;fe,/,m=±l ^ V ' 

+ [A Nm nP{t)Al k - - [Al k A Nmn , p(t) } ) 

N>l;k,m,n=±l ^ ' 

+ ^2 (A Nmn p{t) A ] N , kl - - ^ N , kl A Nmn , p{t) 

N,N'>l:k,l,m,n=±l ^ 
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+ 7n(a; ) 



amy P\t) 



E ( A Lp(t)A ok - ~ [A ok Al 

m=±l ^ 

+ J2 [Al m p(t)A Nlkl -^A N/kl Al 

N'>l;k,l,m=±l ^ 

E [ A N m nP(t)Aok - \ {A kd Nmn , P {t)} 



+ 

N>l;k,m,n=±l 

+ E (^!v m nPW^/-^{^'H^jv m „,PW}) -(23) 

AT,JV'>l;fe,i,m,n=±l ^ ' _ 

The next step will be a comparison between the phenomenological model in the 
form of Eq. (1231) and some suitable forms of the microscopic model in Eq. f fT9|) . 

From a first comparison between Eq. ( jT9l) and Eq. ( 1231) . it is very easy to see that 
the two models coincide when the temperature of the reservoir is zero and its spectrum 
is flat, a condition strictly related to Markov approximation, and under the hypothesis 
that all the imaginary parts of the T coefficients, i.e. all the Lamb shifts, are negligible. 
Indeed under these assumptions all the T coefficients in Eq. (fT9l) are equal to the zero 
temperature spectral density J (see Eq. ([ZED) an d so is half the rate of cavity losses 
7/2 in Eq. (|23|) . Calling 7/2 the common value of the T coefficients, Eq. (jT9l) can be 
rearranged in order to be exactly identical to Eq. ( 1231) with n{ujo) = 0. 

This result provides a condition for the validity of the phenomenological model at 
zero temperature: it is enough that the spectral density is flat over the Bohr frequencies 
involved and that all the renormalization terms arising from the imaginary parts of the 
T coefficients, see Eq. ( 120]) . are negligible. The latter condition on the imaginary parts 
is very often used in the literature, since it is assumed that the Cauchy principal part 
of an integral of the spectral density is very small compared to the real parts of the 
coefficients [16J. The assumption of flat spectrum, instead, depends on the environment 
with which the atom-cavity system is interacting, and on the initial state of the system: 
indeed, since in view of the condition T = the energy of the system cannot increase 
in time, one is allowed to truncate the state space to a subspace with a finite number 
of excitations equal to the maximum number of excitations in the initial state, and it is 
in this space that all the decay rates must be the same. 

Let us examine Eq. (|T9l) from another point of view: the equation is written as 
a sum of many terms, most of which are counter-rotating. Essentially the counter- 
rotating terms are all the terms in the sums wherein the downward jump operator is 
different from the hermitian conjugate of the upward jump operator in the same term 
(for example the term Ao + p(t)Al_ in the first row). What we have previously shown 
with our analysis is that, under the assumption of zero temperature and flat spectrum, 
Eq. (|T9l) can be summed up and can be put in the form of Eq. ([5]) . 

Therefore an aspect emerging from our analysis is that the zero temperature 
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phenomenological model in Eq. (JS]) arises microscopically from an interplay of many 
counter-rotating terms. This point gives an interesting way of looking at the shift in the 
frequency of the Rabi oscillations in the atom-cavity system predicted by Eq. (jSJ) [7j. 
Indeed in the phenomenological model, because of cavity losses, the Rabi frequency is 
not exactly equal to 2fi, but it is equal to 2$7yl — (7/4f2) 2 , i.e., the Rabi frequency 
is shifted by a quantity which is second order in the ratio between the decay rate 7 
and the Rabi splitting Q. The shift may be interpreted as a cooperative effect of the 
counter-rotating terms in the dissipator, analogous to the Bloch-Siegert shift appearing 
in the JC model when counter-rotating terms are added to the atom-cavity interaction 
Hamiltonian |17j . 

When the temperature of the reservoir is different from zero, the comparison 
becomes more difficult to carry on. Let us assume that the imaginary parts of the 
r coefficients are all negligible and that the zero temperature part Jq(u) of the spectral 
density of the environment is flat. According to Eqs. (I20I) and (I2TI) . the real parts of 
all the r(a>) factors, which are proportional to J(uj), contain the quantities n{uS) + 1 
for u > and n(u>) for uj < 0. Therefore in principle one should claim (1) that neither 
can the V coefficients be all equal if T 7^ 0, nor can the V coefficients, and (2) that the 
phenomenological model in Eq. can never be microscopically justified. 

Moreover in this case the microscopic Quasi-RWA master equation may not be cast 
in the Lindblad form. This leads to some problems in its use. On the one hand it may be 
very difficult to solve, since in general there are not any standard techniques developed 
for finding solutions of non-Lindblad master equations. On the other hand there is the 
more fundamental problem that the general form of Eq. (|T9j) may lead to unphysical 
results, since it may violate complete positivity [8]. Hence, what we want to look for is 
a condition under which the phenomenological model is at least approximately valid, so 
that the problem of the cavity losses can be treated by means of Eq. (j3J). 

If during its time evolution the atom-cavity system has at most N excitations, the 
largest Bohr frequency involved in the dynamics is lu + (y/N + 1 + \/~N)£l while the 
smallest one is lu — (y/N + 1 + y/~N)n. Therefore we have to check the magnitude of 
the difference between the rates corresponding to the lowest and to the largest Bohr 
frequency respectively. Since we are assuming that the quantity Jq{uj) in Eq. (j2TI) is 
flat, the difference in the rates is all due to the difference in the photon populations 
n(uj). 

Figure [H shows a plot of the difference in the populations n(u) = 
1/ [exp (tvjj/k-sT) — 1] of two modes of the reservoir with frequencies uj and to + 5uj 
respectively, as a function of the temperature T and of the frequency difference 6~uj. 
From the plot one can argue that, for low temperatures and not too large values of 5uj, 
the error one makes assuming that all the decay rates in the master equation are equal 
(and assuming the same for the excitation rates) can be considered negligible. Since 
the Suj's in Eq. f|T9|) are proportional to V~N, the negligibility of the difference in the 
rates imposes that the number of excitations in the atom-cavity system cannot be too 
large. Therefore we can conclude that, if the difference in the rates can be neglected, 
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Figure 1. Difference An in the average number of photons between a mode of 
frequency lo and a mode of frequency lo + Slo. The frequency lo is held fixed and 
the difference Art is shown as a function of the temperature T (in units of and 
of the frequency difference Slo (in units of lo). 

the phenomenological model in Eq. (J3J) can be considered as a valid model, provided the 
temperature and the number of total excitations are not too large. 

5. Discussion and Conclusive Remarks 

The analysis we have presented shows that, though the microscopic and phenomenolog- 
ical models coincide when the spectrum of the environment is flat, these two models are 
deeply different. The discrepancy, which involves the predictions when the spectrum is 
not flat and/or at T ^ 0, is conceptual anyway, and relies on the fact that the relevant 
physical mechanisms are different. 

The point of view of our approach is the following: the general formalism of master 
equations claims that all the decay and decoherence channels involve the atom-cavity 
system as a whole, and even in the cases wherein it is the cavity only which directly 
decays also for the microscopic model, this may be seen as a cooperative effect of many 
non-resonant terms describing the coupling between populations and coherences relative 
to atom-cavity dressed states. 

A last point is worth mentioning. From our analysis we have understood that 
the success of the phenomenological model relies on the fact that the Quasi-RWA 
microscopic model reduces to it when the spectrum of the environment is flat. The 
situation is quite different when the spectrum is substantially non-flat. This makes us 
interested in the possibility of studying the microscopic model in the latter case, for 
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which a non-Markovian theory may be necessary. This is the scope of our future work. 
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